GLOBAL EXISTENCE FOR SYSTEMS OF NONLINEAR 
WAVE EQUATIONS IN 3D WITH MULTIPLE SPEEDS 



THOMAS C. SIDERIS AND SHU-YI TU 

Abstract. Global smooth solutions to the initial value problem for 
systems of nonlinear wave equations with multiple propagation speeds 
will be constructed in the case of small initial data and nonlinearities 
satisfying the null condition. 



1. Introduction 

This paper is concerned with the Cauchy problem for coupled systems of 
quasilinear wave equations in three space dimensions of the form 

d t V - 4 Au k = C 3 a k p (du)d a dpu j , k = 1, . . . , m, 

subject to suitably small initial conditions. We assume that the propaga- 
tion speeds are distinct, and we refer to this situation as the nonrelativistic 
case. Here, du stands for the full space-time gradient, and CjL(£) = 0(|£|) 
are smooth functions near the origin in M 4m . We shall construct a unique 
global classical solution, provided that the coefficients of the nonlinear terms 
satisfy the null condition. This nonrelativistic system serves as a simplified 
model for wave propagation problems with different speeds, such as nonlin- 
ear elasticity, charged plasmas, and magneto-hydrodynamics. 

The main difficulty in the nonrelativistic case is that the smaller symmetry 
group of the linear operator weakens the form of the invariant Klainerman 
inequality, see Section ^. In order to obtain a viable L°° — 1? estimate 
for solutions, we utilize an additional set of weighted 1? estimates, as has 



been developed in |14], |pC[| , UK- The advantage of this method is the total 
avoidance of direct estimation of the fundamental solution for the linear 
problem or any type of asymptotic constructions. We treat nondivergence 
form nonlinearities which may contain both spatial and temporal derivatives. 

In the 3D relativistic (scalar) case, the null condition was first identified 
and shown to lead to global existence of small solutions in , . Without 
it, small solutions remain smooth "almost globally" but arbitrarily 



small initial conditions can develop singularities in finite time [OJJ19]. Small 



solution always exist globally in higher dimensions [1C], [[Tgl , [pTlJ ] . The 2D 
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relativistic case is rather more complicated. The sharpest results are given 
in p|, but other work appeared previously in M, M. 



The case of nonrelativistic systems has been considered in 3D [22] and 
in 2D Q, Q. We mention also the early work 0, @, |l7| which deals 
with nonresonant interactions. The common theme in these works is the 
direct estimation of the fundamental solution which, as mentioned above, is 
avoided here. 

The statement of the main result is given in section |3| after the introduc- 
tion of some standard notation. The rest of the paper presents the proof. 
To simplify the exposition, we truncate the nonlinearity at the quadratic 
level, but this entails no loss of generality since the higher-order terms do 
not affect the global behavior of small solutions, |11]. 



2. Notation 

Points in R 4 will be denoted by X = (x ,x ,x 2 ,x 3 ) = (t,x). Partial 
derivatives will be written as 8^ = d/dx k , k = 0, . . . ,3, with the abbrevi- 
ations 8 = (80,81,82,83) = (c\, V). The angular-momentum operators are 
defined as 

n = (^1,^2,^3) = x A V, 

where A denotes the usual vector cross product in M 3 , and the scaling oper- 
ator is defined by 

(2.1) S = tdt + rd r = x a 8 a . 

The collection of these seven vector fields will be labeled as 

r = (r 0} ... ,r 7 ) = (d,n,s). 

Instead of the usual multi-index notation, we will write a = (a±, . . . , a K ) for 
a sequence of indices a* € {0, . . . ,7} of length \a\ = k, and 

r a — r . . . r 

Suppose that b and c are disjoint subsequences of a. Then we will say 
b + c = a, if |6| + |c| = \a\, and b + c < a, if |6| + |c| < \a\. 
The d'Alembertian will be used to denote the operator 

□ = Diag(D 1 , ... ,n m ) with D k = 8 2 t - c\ A. 

For convenience, we will assume that the speeds are distinct 

Ci > • • • > Cm > 0. 

It is also possible to treat the case where some of the speeds are the same, 
see the remark following the statement of Theorem 3.1. This operator acts 
on vector functions u : M 4 — > M m . The standard energy then is defined as 



m 

E l (u(t)) = J2 I [\d t u k (t,x)\ 2 + 4\Vu k (t,x)\ 2 }dx, 

k=l ■ 
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and higher order derivatives will be estimated through 
(2.2a) E K (u(t))= Ei(F a u(t)), « = 2,3,... 

|a|<K-l 

In order to describe the solution space, we introduce the time-independent 
vector fields A = (Ai, . . . , A7) = (V, Q, rd r ). Define 

H£(M. 3 ) = {/ G L 2 (M 3 ;M m ) : A a f G L 2 , \a\ < «}, 

with the norm 

(2.2b) \\f\\ H «= IIAVIIl- 

\a\<K 

Solutions will be constructed in the space Hfi(T) obtained by closing the 
set C oo ([0,T);C^(R 3 ,R m )) in the norm sup E l J 2 {u(t)). Thus, 

0<i<T 

{K-l 
u(t,x):du(t,-)e fl^aCT);^- 1 ^ 
3=0 



By flOD , it will follow that ff«(T) C C K - 2 ([0,T) x M 3 ;M m ). 

An important intermediate role will played by the weighted norm 

m 

(2.2c) X K (u(t))=^Yl E \\(ckt-\x\)d a T b u k (t)\\ L , m , 

k=l \a\=2 \b\<K-2 

where we use the notation (p) = (1 + Ipl 2 ) 1 / 2 . 

3. Main Result 

Consider the initial value problem for a coupled nonlinear system of the 
form 

(3.1) \3u = N(u,u) 

in which the components of the quadratic nonlinearity depend on the form 
(3.2a) N k (u, v) = C^d^dpd^ . 

Summation is performed over repeated indices regardless of their position, 
up or down. Greek indices range from to 3 and Latin indices from 1 to m. 

Existence of solutions depends on the energy method which requires the 
system to be symmetric: 

The key assumption necessary for global existence is the following null 
condition which says that the self-interaction of each wave family is nonres- 
onant: 

(3.2c) C k f 1 X a X p X 1 = for all Ie4 k = 1, . . . , m, 
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with the null cones 

o4 . v2 



Nfc = {X € R : - + ^ + X 3 ^) = 0}. 

Theorem 3.1. Assume that the nonlinear terms in (3.2a) satisfy the sym- 



metry and null conditions (3.2b), ( p. 2c ) . Then the initial value problem for 



(|3.l[) with initial data 



d a u(0) € iZ£ _1 (R 3 ), k>9 



satisfying 

(3.3) E l J 2 2 (u{Q)) exp CE^(u(0)) < s, 

with e sufficiently small, has a unique global solution u G Hfi(T) for every 
T > 0. The solution satisfies the bounds 

E l J^ 2 (u{t)) < 2s and E K (u(t)) < AE K (u(0))(tf £ . 

Remark. We briefly discuss the case when some of the speeds are repeated. 
Suppose that only £ < m of the speeds, c\ > Ck 2 > ■ ■ ■ > Ck e are distinct. 
For p = 1, . . . ,£, let I p = {k : 1 < k < m, ct = Cfc p }. The null condition is 
now extended to be 

&* 1 X a XpX 1 = 0, for all Ie\, (i, j, k) G I 3 , p = 1, . . . , L 

The proof can easily be adjusted to handle this more general case. 

4. Commutation and Null Forms 

In preparation for the energy estimates, we need to consider the commu- 
tation properties of the vector fields V with respect to the nonlinear terms. 
It is necessary to verify that the null structure is preserved upon differenti- 
ation. 



Lemma 4.1. Let u be solution u of ( p.ip in Hft(T). Assume that the null 
condition ( |3.2cj ) holds for the nonlinearity in (3.2a). Then for \a\ < k — 1, 



nT a u= N d (T b u,T c u), 

b+c J rd=a 



in which each Nd is a quadratic nonlinearity of the form (3.2a) satisfying 



(3.2c) . Moreover, if b + c = a, then N d = N . 
Proof. First we note the well-known facts that 

[d, □] = 0, [Q, □] = 0, [S, □] = -2D. 
Recalling the definition ( p.2a[) , we set 

[T, N](u, v) = TN(u, v) - NiTu, v) - N(u, Tv). 



This is a quadratic nonlinearity of the form ( 3.2a ). Thus, if [r,iV] is null 
for each T, then the result follows by induction. In fact, if d = (di, . . . , dk), 
then Nd is the fc-fold commutator Nd = [Td k , [■■■■, [T^, N]]]. 
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A simple calculation shows that 

[d, N] (u, v) = and [S, N] = -3N{u, v) . 

Thus, these commutators are null if N is null. 

We can express the angular momentum operators as Q\ = ex^x^dy, 
A = 1,2,3, where £\^u is the tensor with value +1, —1 if \\iv is an even, 
respectively odd, permutation of 123, and with value otherwise. Using 
this, we find that the kft 1 component of N] is 

[Q x ,N\ k (u,v) = d i ^ r d a ^dp& r v k 

with 

To see that this commutator is also null, write 

h k (X) = C kk ^ k X a XpX 1 and h k (X) = C k p k X a XpX 1 . 

Then h k (X) = Dh k (X)Y x with = e x ^X u . Now the null condition says 
that h k (X) = for X £ N^. But since Y x is tangent to at X, we have 
h k {X) = for X e N fc . This implies that [fi A , N] is null. □ 



5. Estimates for Null Forms 

The utility of the null condition is captured in the next lemma. The 
presence of the the terms with the weight (c^t — r) in these inequalities is 
explained by the absence of the Lorentz rotations in our list of vector fields 

r. 

Lemma 5.1. Suppose that the nonlinear form N(u,v) defined ( |3.2a| ) satis- 
fies the null condition (3.2c). Set cq = min{cfc/2 : k = 1, ... ,m}. For u, v, 
w G C 2 ([0,T] x M 3 ;lR m ) and r > c t, we have at any point X = (t,x) 

(5.1a) \c k a kk d a u k d p d lV k \ 



< 



c r 



(X) 



|lV||dV| + \du k \\dTv k \ + (c k t - r)\du k \\d 2 v ki 



and 



(5.1b) \C$frd a u%v%w k \ 



< 



C 
(X) 



\r u k \\dv k \\dw k \ + \du k \\rv k \\dw k \ + \du k \\dv k \\r w k \ 

+ (c k t-r)\du k \\dv k \\dw k \ 



which (X) = (1 + \X\ 2 fl 2 . 



6 THOMAS C. SIDERIS AND SHU-YI TU 

Proof. Spatial derivatives have the decomposition 

x 



v = -d r 

r 



An. 



So if we introduce the two operators D k = | (d t ± c k d r ) and the null vectors 
Y k = (l,±x/c k r) £ 3Sffc we obtain 



(5.2) av) = (p fe -+w 

On the other hand, if we write 



0, 



Afl . 



c k r< 



K c k t + r c k t + r K 
the formula fl5.2|) can be transformed into 



Cht — r , 



c k t + r 



c k t + r 



c k r^ 



Thus, we have 
(5.3a) 



9=y k -D k +R 



Now, we may assume that \X\ > 1, for otherwise the estimates are trivial. 
But then it follows that 1/r and \/{c k t + r) are bounded by C/(X), and as 
a consequence we have 

(5.3b) \Ru\ < C(X) -1 [|ru| + (c k t - r)\du\] 

Using ( |5.3a ), we have 



(5.4) CSfrdavfdfoJ = C^[Y k - a Y k - p Y k - lD 7u\D~ k fv k 

+ R a u%d lV k + Y^D-^Rpdy + Y h -D-u k Y- R D-R y v k }. 



ka^k A k/3^k 



The first term in (5.4) vanishes since N obeys the null condition, and by 
( |5.3b ) the remaining terms in ( |5.4[) have the estimate ( [5. la ). 
The proof of ( 5.1b ) is similar. 

□ 



6. Sobolev Inequalities 

The following Sobolev inequalities involve only the angular momentum 
operators since we are in the nonrelativistic case. The weight {ct — r) com- 
pensates for this. We use the notation defined in ( [T2p , ( |2~2b| ), (|]§). 
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Lemma 6.1. Let u G H£(T), with 1 K {u{t)) < oo. 

(6.1) {r)V 2 \T*u(t,x)\<CEV*(u(t)), 

(6.2) (r)\dT a u(t,x)\ <CEl/ 2 {u{t)), 

(6.3) {r){ct-r) l / 2 \dT a u{t,x)\ 



\a\ + 2 < k 
\a\ + 3 < k 



(6.4) 



< C El ,2 {u{t)) + % R (u{t)) , |a| + 3 < k 
(r){ct-r)\d 2 r a u{t,x)\ < CT K {u{t)), \a\ +4 < k. 
Proof. This result is essentially Proposition 3.3 in [21|. 

7. Weighted Decay Estimates 



□ 



The main extra step in the nonrelativistic case is to control the weighted 
norm X K (u(t)). This will be accomplished in this section by a type of boot- 
strap argument. 

Lemma 7.1. Let u G H£(T). T/ten 

(7.i) x>(t))<c £y 2 Kt))+ E ll(i + Onr a ^)ll| 2 

|a|<K-2 

Proof. Recall that the weighted norm involves derivatives in the form d 2 T a u 
In the case when d 2 = V9, the result was given in Lemma 3.1 of 0. 
Otherwise, if d 2 = d 2 , then the result is an immediate consequence of (2.10) 
in 0. 

Now we assume that u solves the nonlinear PDE. 

Lemma 7.2. Let u G Hfi(T) be a solution of ([O]). Define k! 
Then for all \a\ < k — 2, 



□ 



[*fl]+3. 



(7.2) \\{t + r)UY a u{t)\\ 2 L2 



< C[X K ,{u{t))E l J 2 {u(t)) + T K {u{t))E l J, 2 {u{t))\- 



Proof. By Lemma |4jj, we must estimate terms of the form 

\\( y t + r)dr b u i d 2 r c ui\\ L 2, 

but since (t + r) < C(r)(cjt — r), we will consider 

(7.3) \\{r)(cjt - r)dT b u l d 2 r c ui\\ L 2, 

with b + c < a, and \a\ < k — 2. 

Let m = [^-^-] = k' — 3. We separate two cases: either \b\ < m or 



|c| < m — 1. In the first case, Q7.3j ) is estimated as follows using ( |6.2| ): 

||(r)arV|| i0 c||(c,i-r)9 2 rV|| L2 < C^ 2 (n(t))X K (n(t)). 
Otherwise, we use (|6.4| ) to estimate ( |7.3| ) by: 

||9rV|| L 2||(r)( Ci t-r)5 2 rV ||ioc < C^ /2 (n(t))X K /(u(t)). 



s 
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□ 

The next result gains control of the weighted norm by the energy. We 
distinguish two different energies, the smaller of which must remain small. 
In the next section, we will allow the larger energy will grow polynomially 
in time. 



Lemma 7.3. Let u € Hfi(T), k > 8, be a solution of ( p.l|) . Define = 
k — 2, and assume that 

e = sup E]J 2 {u{t)) 

0<t<T 

is sufficiently small. Then for < t < T, 

(7.4a) X M (u(i)) < CE]/\u{t)) 

and 

(7.4b) 1 K iu(t)) < CE^\u(t)). 

Proof. Let n' = + 3, /i = n — 2. Since fi > 6, we have // < /x. Thus, 

by Lemmas 7.1 and [7.2| , we find using our assumption 

X„(u(*)) < C[E l J\u{t)) +e X„(u(t))]. 



Thus, if £q is small enough, the bound ( 7.4a[) results 



Again since k > 8, we have k' = ] + 3</i = K — 2. From Lemmas 



7. 1 and |7.2j we now have 

X„(«(t)) < C[Ell 2 (u(t)) + X^u^El' 2 ^)) + X K {u(t))Ey 2 {u{t))]- 
If we apply ( 7.4a| ) and our assumption, then 

X K (n(t)) < C[El/ 2 (u(t)) + e X K (u(t))}, 



from which ( |7.4b ) follows. □ 



8. Energy Estimates 
General energy method. In this section we shall complete the proof of 



Theorem 3J. Assume that u(t) G Hfi(T) is a local solution of the initial 
value problem for (|3.l| ). Our task will be to show that E K (u(t)) remains 
finite for all t > 0. To do so, we will derive a pair of coupled differential 
inequalities for (modifications of) E K (u(t)) and £^(u(t)), with fj, = k — 2. If 

I K 1/2 

(I3J) holds then (u(0)) < £. Suppose that Tq is the largest time such 



1/2 

that Efi (u(t)) < 2s, for < t < T with e small enough so that Lemma |7.3| 
is valid. All of the following computations will be valid on this time interval. 
Following the energy method, we have for any v = 1, . . . , k, 

KHt))= Y, [(aT a u (t),d t r a u(t))dx, 

\a\<v-l 
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.1) K(u{t))= Y, E f {N d {T b u,T c u),d t T a u)dx. 



\a\<u— 1 b+c+d=a 



Terms in (|S,lD with b = 0, c = a, and |a| = v — \ are handled with the 



aid of the symmetry condition ( 3.2b ) which allows us to integrate by parts 
as follows. Recall that from Lemma 4.1, Nd = N when b + c = a. 



(N(u,T a u),d t T a u)dx =&H \ dau i d/3d 1 u j d t u k dx 



= C a^j d y [d a u i d fi r a v>d t r a u k \dx 
~ Cft / dad^dpT^dtT^dx 
~ C T Pl j d a u%T a yPd t d y T a u k dx 

= C fpo d t j d a u l dpT a uid t T a u k dx 
~ C Tm j d a d 7 u%T a uid t r a u k dx 
- \ C Sj[ d a u%[d fl r a v?dL r r a u k ]dx 

=\ ^VjS dt J d a u l d p T a uid s T a u k dx 
~ C Tm j d a d^d p T a u^d t T a U k dx 

+ \ j dtd«U%T a V>d i r a U ]l dx, 

using the symbol ry 7< 5 = Diag[l, — 1, — 1, — 1]. The first term above can be 
absorbed into the energy as a lower order perturbation. Define 



E v (u{t)) = E u (u{t)) ~\ Y C 2W / d a u%T a u j d s T a u k dx. 

\a\=u—l 



The perturbation is bounded by C\\Vu\\L°°E v (u(t)), but by (|6.2[) , the max- 

1/2 1/2 

imum norm ||Vit||i°o is controlled by E 3 < E^ < 2e. Thus, for small 
solutions we have 



J.2 



(l/2)E u (u(t)) < E u (u(t)) < 2E u (u(t)). 
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Returning to fl3,l| ), we have derived the energy identity 
(8.3) K(u(t))= E f {N d {T b u,T c u),d t T a u)dx 

\a\<v-l 6 l + l c + d = a 



\a\=u—l 



J (N(T b u,r c u),d t T a u)dx 

b+c=a 

~ Cf^ J d a d 7 u%T a ^d t T a u k dx 
+ \ C t^ J dtd^dpT^d^dx 



Higher energy. For the first series of estimates we take v = k in 
We obtain immediately 

(8.4) E' K {u{t))<C E E W b u l d 2 T c ui\\ L 4dT a u k \\ L ,. 



SJ). 



i,j,k \a\<K— 1 



b+c<a 



In some cases, the indices i and j have been interchanged. In the sum on 



the right-hand side of (^4), we have either \b\ < k' or |c| < «' — 1, with 
= [?|] . Note that since k > 9, we have k' + 3<k — 2 = /i. We will also 

use that (t) < C(r)(cjt - r). 

In the first case, we estimate using (|6.2| ) and 7.4b| 



l|sr 6 « i 5 2 r c u 3 " || L2 ^c^r 1 1| (r)dr b u i \\ Lao \\ (c,-t - r)s 2 rv \\ L2 

<C(i)- 1 S | 1 £ 3 (u(t))X (t («(t)) 
In the second case, we use ( |6.4| ) and then ( [T.4aj) 

I [ ar V a 2 rv 1 1 L2 < c (t } - 1 1 1 ar V 1 1 L2 \ \ (r } ( Cj t - r ) d 2 r c u j \ \ L °° 
<c(t)- 1 J ey 2 ( u (t))X| C | +4 (n(t)) 

<C(t)^Ell\u{t))X»{u{t)) 
<C{t)- x El'\u{t))Ey\u{t)). 

Going back to ( |8.4| ) and recalling ( |8.2[ ) , we have established the inequality 
5.5) # B («(t)) < C^Ey^u^EMt)) 
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Lower energy. The second series of energy estimates will exploit the null 
condition. We return to fl3,3| ) now with v = \i = k — 2. The resulting 
integrals on the right-hand side of ( |8.3|) will be subdivided into separate 
integrals over the regions r < cot and r > c^t. Recall that the constant cq 
was defined in Lemma |5.1| . 



Inside the cones. On the region r < c^t, we have that the right-hand side of 
( |8.3[) is bounded above by 



Y Y Y l|5rV5 2 rV5r a u fc || L i 



(r<c t)- 



i ,j,k b+c<a b+c<a 



Since r < c^t, we have that [cit — r) < C(t) for each i = 1, 
using (ETH), a typical term can be estimated by 



, m. Thus, 



C(t)-V 2 \\{cit - r^dT^icjt - r)d 2 T c u j dT a i 



\L 1 (r<cot) 



< C(t)-^' 2 \\{cit - r) ir2 dT b u l \\ L oo\\( Cj t - r)d 2 r c u j \\ L 2\\d t r a u k \\ L 2 



\- 3 /2||/,, / _^l/2 ar 6..<, 

< C(t)- 3 / 2 [^ 3 (n(t)) + X| fe | +3 (u(t))]X| c | +2 (n(t))4/ 2 (n(t)). 

In the preceding, we have \b\ + 3 < k, \c\ + 2 < /i, and |a| + 1 < //. With the 
aid of Lemma 7.3, we have achieved an upper bound of the form 

c<r 3/2 W0)^ /2 («(*)) 

for the portion of the integrals over r < c$t on the right of ( |8.3[ ). 



Away from the origin. It remains to estimate the right-hand side of 
for r > cot. 

First, we consider the nonresonant terms, i.e. those for which (i,j,k) ^ 
(k, k,k). If i ^ j and r > c t, then (t) 3 / 2 < C{r)(at - r) 1 / 2 (c i t - r). Using 
( |6.3|) we have the estimate 

l|srWrvarV|| il(r > cot) 

< C(t)- 3 / 2 \\(r)( Ci t - r) 1 ' 2 dY b u i \\ L o a \\{c j t - r)d 2 T c u j \\ L 2\\dT a u k \\ L 2 
< C{ty 3 ' 2 [E^ 2 +3 (u(t)) + X w+3 (u(t))] l lcl+2 (u(t))E^ +1 (u(t)) 

<C(t)- 3 / 2 ^(n(t))i?y 2 (n(t)). 

Otherwise, if j ^ k, we pair the weig ht (r)(c fc t-r) 1 / 2 with <9r 2 u k in L°° to 
get the same upper bound. 

We are left to consider the resonant terms in Q8.3| ), i.e. (i, j, k) = (k, k, k), 
in the region r > c$t. It is here, finally, where the null condition enters. An 
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application of Lemma 5.1 yields the following upper bound for these terms: 

c^y'H E [l|r b+ Wr^arV|| L1{r > cot) 



k b+c=a 



+ \\dT b u k dT c+l u k dT a u k \\ L , {r > Cot) 



+ \\( Ck t - r)dr b u k d 2 r c u k dr a u k \\ LHr > C0t) 



We still need to squeeze out an additional decay factor of (t)' 



-1/2 



Since r > cot, we have (r) < C{t). Thus, we have using (6.1) 



\\T b+1 u k d 2 T c u k dT 2 u k \\ LHr > C0t) 

< C{ty^E]l^{u{t))E^u{t)) 

< C{t)- l l 2 El/\u{t))EMt))- 

In a similar fashion, the second term is handled using ( |6,2[ ): 

\\dT h u k dT c+l u k dT 2 u k \\ LHr > C0t) 

< C{t)- l E\^{u(t))E,{u(t)) 
<C{ty l E l J 2 (u(t))E,(u{t)). 



The final set of terms are estimated using ( p.2[) again and ( 7.4a ). 

\\( Ck t - r}dr b u k d 2 r c u k dr a u k \\ LHr > C0t) 



Combining all the estimates in this subsection, we obtain, thanks to (8.2), 
the following inequality for the lower energy: 



£>(i)) < C(t)-*l 2 E,{u{t))Ell 2 {u{t)) 
<C{t)-^ 2 E,(u{t))Ell 2 {u(t)). 



(8.6) 



Conclusion of the proof. By (|8.2|), we have that the modified energy 



"1/2 

satisfies Eft («(*)) < Ce for < t < Tq. So from ( |8.5| ) , we find that 

E K (u(t)) < E K (u(0)){t) C£ , 
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provided e is small. Inserting this bound into Q3,6| ) and using ( |8.2| ), we 
obtain 



(l/2)E a (u(t)) < E^{u{t)) < E^umexpCIE^iuiO)) 

< 2E^u{t))exp2CIE 1 J 2 (u(0)) < 2s 2 , 



with / = /q°°(s) 3 / 2+Cs ds. With this we see (u(t)) remains strictly lees 
than 2e throughout the closed interval < t < Tq. This shows that E K (u(t)) 



is bounded for all time, which completes the proof of Theorem 3.1 
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